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ABSTRACT 

We reexamine the age distribution of star clusters in the Antennae in the context of N- 
body+hydrodynamical simulations of these interacting galaxies. All of the simulations that account 
for the observed morphology and other properties of the Antennae have star formation rates that vary 
relatively slowly with time, by factors of only 1.3 — 2.5 in the past 10 8 yr. In contrast, the observed 
age distribution of the clusters declines approximately as a power law, dN/dr oc t 7 with 7 = —1.0, 
for ages 10 6 yr < r < 10 9 yr. These two facts can only be reconciled if the clusters are disrupted 
progressively for at least ~ 10 8 yr and possibly ~ 10 9 yr. When we combine the simulated formation 
rates with a power- law model, / SU rv oc r s , for the fraction of clusters that survive to each age t, we 
match the observed age distribution with exponents in the range —0.9 < 5 < —0.6 (with a slightly 
different S for each simulation). The similarity between S and 7 indicates that dN/dr is shaped mainly 
by the disruption of clusters rather than variations in their formation rate. Thus, the situation in the 
interacting Antennae resembles that in relatively quiescent galaxies such as the Milky Way and the 
Magellanic Clouds. 

Subject headings: galaxies: individual (NGC 4038/39) — galaxies: interactions — galaxies: star 
clusters: general — methods: numerical 



1. INTRODUCTION 

Interacting galaxies in the nearby universe are labo- 
ratories for direct studies of several physical processes 
that were important in the formation and early evolu- 
tion of galaxies. From such studies, we hope to learn, 
for example, how interactions and mergers affect the cy- 
cle in which baryonic matter is converted from diffuse 
interstellar gas into dense molecular clouds, then into 
star clusters, and eventually, by disruption, into a rel- 
atively smooth stellar field. It is clear that interactions 
and mergers boost the rate of star and cluster formation. 
But do they also change the rate at which clusters are 
disrupted? This is the question we address in this paper. 

The most intensively studied interacting galaxies are 
the Ant ennae (NGC 4038/3 9), at a distance D ~ 
20Mpc ([Schweizer et al.ll2008| ). They consist of two nor- 
mal disk galaxies that began to collide a few x 10 8 yr 
ago. The stellar population and interstellar medium 
(ISM) of the Antennae have been observed over an enor- 
mo us range of wavelengths, from X-ray to radio (see , 
e.g.lZhang et al.ll200lHHibbard et al.l | 2001t jKassin et al 



2003: IZezas et al.| ]2006: Br andl et all 120091 : iKlaas et al 
20101 and references therein). The star clusters have 



been the focus of numerous studies based on observa- 
tions with the Hubble Space Telescope (HST), culminat- 
ing in well-deter mined luminosity, mass, age, and space 
distributions (see IWhitmore et al.l I2010L and references 
therein) . 

The Antennae have also be en the focus of seve r al dy- 
namical simulations, first by iToomre &: Tpomrd ([1972L 
hereafter TT72) and then bv lBarnesI (|1988h . These pio- 
neering studies demonstrated that gravity alone can ac- 
count for the gross features of the observed morphol- 
ogy and kinematics of the stellar components of the 



merger. Subsequent simulations have included an inter- 
stellar medium and star formation, with the additional 
goal of matching the ob served space distribution of young 
stars in the Antennae ( jMihos et al.lll993l : iTeyssier et al.l 
[20101 : iKarl et"aTl MHI hereafter MBR93, TCB10, and 
K10, respectively). There have also been two recent at- 
tempts to match the observed age distribution of the 
clusters, with different assumptio ns about their disrup- 
tion histories ([Bastian et al.ll2009l K10). 

The purpose of this paper is to reexamine the is- 
sues raised by the observed age distribution of the clus- 
ters in the Antennae. In Section [2j we review the evi- 
dence for a quasi-universal age distribution of star clus- 
ters in different galaxies, and in Section [3j we assem- 
ble the star formation histories from all available N- 
body+hydrodynamical simulations of the Antennae. We 
then combine these and compare the results with ob- 
servations in Section [4[ We summarize our conclusions 
and their implications in Section [5] In particular, we 
show in this paper that there is nothing special about 
the disruption history of clusters in the interacting An- 
tennae galaxies; it is similar to that in quiescent (non- 
interacting) galaxies. 

Before proceeding, we offer a few remarks about 
nomenclature. We use the term "cluster" for any con- 
centrated aggregate of stars with a density much higher 
than that of the surrounding stellar field, whether or 
not it also contains gas and whether or not it is grav- 
itationally bound. This is the stand ard definition in the 
star formation communi ty (see, e.g. lLada fc Ladall2003l : 
iMcKee fc Ostr iker 2007). Some authors use the term 
"cluster" only for gas-free or gravitationally bound ob- 
jects. Such definitions are not appropriate in the present 
context for two reasons: (1) A key element in our analy- 
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sis is the connection between the formation rates of stars 
and clusters. We would break this connection artificially 
if we were to exclude the gas-rich clusters in which stars 
form. (2) It is virtually impossible to tell from observa- 
tions which clusters satisfy the virial theorem precisely 
and which do not, especially at the distance of the An- 
tennae. Indeed, iV-body simulations show that an un- 
bound cluster retains the appearance of a bound cluster 
for remarkably long times, m ore than 10 crossing times 
(jBaumgardt fc Kroupall2007l) . 

2. DISRUPTION OF STAR CLUSTERS 

Star clu sters form in the dense inner parts of molecu - 
lar clouds (|Lada fc Ladall2003HMcKee fc Ostrikerll2007D . 
Most of them are subsequently destroyed by various 
mechanisms, beginning with the expulsion of interstel- 
lar material by massive young stars ("feedback"), later 
mass loss from intermediate- and low-mass stars, tidal 
disturbances from passing molecular clouds, an d stellar 
escape driven by internal two- body relaxation (|Spitzerl 
119871: iBinnev fc Tremainell2008l ). This leads to the even- 
tual dispersal of stars from the clusters into the sur- 
rounding stellar field. In the Milky Way, the fraction 
of stars in clusters declines with age r, from / c i us > 
50% at r ~ 10 6 yr to f r i,, s < 1% at r ~ 10 9 yr 
(|Binnev fc Merrifieldl [19981 ILada fc Ladal 12001 . The 
goal of this paper is to test whether a similar situation 
holds in the Antennae. 

The age distribution of clusters dN/dr in a galaxy rep- 
resents the formation rate (dN/dr)[ mm modified by sub- 
sequent disruption, leaving behind the survival fraction 
/survM of clusters at each age r: 



dN/dr = f sut y (dN/dr) Ion 



(i) 



Thus, if the formation rate varies slowly with time, the 
age distribution primarily reflects the disruption history 
of the clusters, i.e., dN/dr ^ /surv If, on the other 
hand, the survival fraction varies slowly, the age dis- 
tribution primarily reflects the formation history, i.e., 
dN/dr^ (dN/dr) f orm - We assume throughout this pa- 
per that the formation rates of clusters and stars, 
(dN/dr){ mm and dN^/dr, track each other: 

(dN/dr)f OTm = c ■ dN„/dr with c = constant. (2) 

This is certainly true if most stars fo r m in clusters, as in 
the M ilky Way (|Lada fc Ladal 120031: iMcKee fc Ostrikerl 
2007). Equation ([2]) is the most direct connection be- 
tween (dN / dr)i ovm and dN^/dr; hence it may be a good 
approximation whether or not most stars form in clus- 
ters. In the Antennae, we know from Ha observations 
that at least 20% and possibly all stars form in clusters 
(iFall et al]l2005h. 

iFall et al.l (|2005l ) found that the age distribution of 
massive clusters in the Antennae galaxies can be ap- 
proximated by a power law, x(r) = dN/dr cx r 7 with 
7 w —1. The mass function of the clusters is also 
a power law, MM) = dN/dM cx with /3 « -2 
([Zhang fc Falll fl999f ) . In fact, these results are parts of 
a broader finding; the bivariate distribution of masses 
and ages can be approximated by a product of power 
laws: g(M,r) cx ip(M)x{r) cx M^r 7 . This model has 
been derived from observations of massive clusters in the 
Antennae (roughly M > 10 4 (t/10V) - 7 M q ). The de- 
composition of g(M, r) into a product of ij)(M) and x( T ) 



implies that the (net) formation and disruption rates 
of the clusters are independent of their masses 1 . The 
power-law model for g(M,r) has been confirmed in sub- 
sequent observational studies of the Antennae clusters 
(iWhitmore et al.l [20071: [Fall et al.l [20091 iWhitmore et al.l 
20 10). 

Similar results have now been obtained for the age 
distributions of clusters in about 20 other galaxies (al- 
though with smaller samples and thus larger uncer- 
tainties than for the Antennae galaxies ) . These in- 
clude the Milky Way (ILada fc Ladal I2003D : the Larg e 
and Small Magellanic Clouds (Chanda r et al.l [2010a); 
NGC 1313, NGC 4 395, NGC 523 6 (M83 ), NGC 7793 
Mora et al.l [200l [Chandar et al.l l20l0ch: NGC 92 2 



Pellerin et al.ll2010l): NG C 3256 dGoddard et al.l l20Tol): 
Arp 284 ([Peterson et al.ll2009D : and nine nearby dwarf 
galaxies (jMelena et al.ll2009ft . The masses and ages of 
the clusters in these galaxies cover the ranges 10 2 M Q < 
M < 10 6 Mq and r < 10 9 yr. The galaxies themselves 
are also diverse, ranging from dwarf to giant, quiescent 
to interacting. Yet in all cases, the observed age dis- 
tribution of the clusters can be represented by a power 
law, 



dN/dr 



with 



l-0<7<-0.7, (3) 



wh ere the exponen t s have uncertainties A7 w 0.1 — 0.2. 
See IChandar et al.l (|2010bD for a more complete discus- 
sion of these results and the methods used to obtain 
them. 

The stellar age distribution is known in detail only 
for t hree of the galaxies menti oned above: the Milky 
Way (Binney & Merrifield 1998, and references therein), 
the LMC ([Harris fc Zaritskvl 120091 ) . and the SMC 
([Harris fc Zaritskvl l2004f) '. In these cases, the star for- 
mation rate (SFR) has been constant to within a factor 
of 2 over the past ~ 10 9 yr, and we can be certain that 
the observed age distribution of the clusters primarily 
reflects their disruption history. From Equations ([TJ and 
^ and the observed dN/dr we obtain 2 



,/surv OC r 



with 



7 



-1. 



(4) 



Less is known about the star formation histories in the 
other galaxies. For some, the rate may have increased 
over the past ~ 10 9 yr, while for others, it may have de- 
creased. However, if we assume that the average forma- 
tion rate for the sample (which includes a fairly typical 
mix of galaxies) has been roughly constant, then Equa- 
tions ([1]) and ([2]) and the observed dN/dr again imply 
/surv 2 t S with S « 7 « — 1. In fact, the average SFR 
over all galaxies has declined slightly i n the past ~ 10 9 yr 
([Lilly et all 119961: iMadau et al.l ll996D . The alternative 
to this interpretation is that all 20 galaxies have the 

1 There is no evidence for mass-dependent disruption of clusters 
in the Anten nae. Claims for m ass-dependent disruption in other 
galaxies (e.g., Giclcs ct al. 2005) have al l been contradicted b y bet- 
ter data and more direct analyses (e.g.. IChandar et al.H201ll) . 

2 Strictly speaking, the survival fraction must have the limiting 
behavior / sur v — > 1 for r — ¥ (i.e., before any clusters are dis- 
rupted), whereas the power-law model diverges. We could fix this 
by introducing a bend at some young age (< 10 6 yr) such that 
/surv ~ 1 for r < tj, and / surv oc t s for r > t^. This is not neces- 
sary in practice, however, because we only make comparisons with 
binned data (in Section 4), for which the exact behavior of / sur v 
near r = is irrelevant. 
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Fig . 1 . — Evolution of the star formation rates in different Af-body+hydrodynamical simulations of the Antennae galaxies. B lue lIKarl et al 



201C , K10, and this paper): strong feedback (dashed), intermediate feedback (solid), weak feedback (dot-dashed); green IITevssier et al 



20101 . TCB10): low-resolution (dotted), high-resolution (dot-dashed); and red dotted llMihos et a ll 1993. MBR93). The time of best match 
for each simulation is t = (at the vertical solid line). The time of first pericenter passage is indicated by colored vertical bars. Straight 



horizontal lines show the assumed constant star formation rates in the pre-interaction disks. 

ested reader is referred to lFall et "all f|2009T ) and lFall et al.l 
(|2010h for further discussion of disruption mechanisms 
and how they relate to the mass and age distributions. 



same rising star formation rate. However, the proba- 
bility of this happening by chance is minuscule (roughly 
(1/2) 20 « 1CT 6 for equal numbers of rising and falling 
SFRs). 

For the Antennae, we have another indication that 
the observed age distribution of the clusters reflects 
their disruption history rather than their formation his- 
tory. Specifically, dN/dr has approximately the same 
power-law shape in r egions separated by ~ 10 kpc 
(Whit more et al.1 12007;. This observation can be ex- 
plained, in principle, either by disruption or by syn- 
chronized formation throughout the Antennae. Indeed, 
the spatial uniformity of dN/dr, combined with causal- 
ity restrictions, places interesting constraints on pos- 
sible temporal variations in the formation rate as fol- 
lows dFall et all I2005LI2003 I Whitmore et all 120071 ) . The 
timescale for variations in the formation rate, if these 
are driven by gravitational interactions, is the galactic 
orbital period, ~ 10 8 yr. Similarly, the communication 
time is 10 8 yr or 10 9 yr for a signal traveling a distance 
of 10 kpc at a velocity of 100 km s _1 or 10 km s _1 , re- 
spectively, plausible values for pressure disturbances in 
the ISM. In either case, we expect the formation rate to 
vary relatively slowly, an expectation borne out nicely 
by the A^-body+hydrodynamical simulations presented 
in the next section. 

The mechanisms that disrupt clusters include the 

following: (1] re mova l of ISM by stellar feedback, 

t < 10 7 yr (Hills 1980); (2) continued mass loss from 
intermediate- and low-mass stars, 10 7 yr < r < 10 8 yr 
(|Chernoff fc Weinberg! Il990[ ): (3) tidal d isturbances b y 
passing molecular clouds, r > 10 8 yr (Spitzer 1958); 
(4) stellar escape driven by i nternal two-body relax- 
ation, t > 10 9 yr (|Spitzerlll987l ). The timescales quoted 
above are highly approximate, and some of the mech- 
anisms must in fact operate simultaneously. The inter- 



In the present context, it is important to note that mech- 
anisms (1), (2), and (4) have little or no dependence on 
the properties of the host galaxy. Only mechanism (3) — 
encounters with molecular clouds — is expected to have 
such a dependence, through the mean density of molec- 
ular gas. In practice, however, this is likely to make rel- 
atively little difference over much of the observed range 
of ages, and we therefore expect / SU rv to be similar in 
different galaxies, consistent with the observations sum- 
marized above. 

3. Af-BODY+HYDRO DYNAMICAL SIMULATIONS 

Interacting and merging galaxies like the Antennae 
generally show signs of enhanced star formation in the 
past few x 10 8 yr, along with more complex morphology 
and kinematics, compared with isolated disk galaxies. 
In this Section, we briefly present and analyze all pub- 
lished and two new simulations of the Antennae galaxies 
that make specific predictions for their star formation 
history. We do not claim that any of these simulations 
uniquely represents the real Antennae system. In par- 
ticular, the spatial and temporal resolution in the simu- 
lations is much too low to model the ISM accurately on 
the scales of individual star clusters dBekki et al.l 120021 : 
iKravtsov fc Gnedinl [20051: |UeiaL|l2004j) ;:i . Nevertheless, 
these simulations provide us with a suite of plausible star 
formation histories for the Antennae. We adopt this ap- 
proach because it would be difficult, if not impossible, 
to determine the stellar age distribution in the Antennae 
directly from observations. 

3 See[Bournaud et al. (2008) fo r a first attempt to s imulate clus- 
ter formation directly and lRenaud ~etaD. (2008, 20Q3) for the pos- 
sible effects of compressive tides on cluster formation. 
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Fig. 2. — Star formation history from K10 with intermediate feed- 
back, g oqs = 0.5 (solid line, see also Figure[T}. The coarse-grained, 
fluctuating age distribution of the spawned stellar particles from 
Figure 4 of K10 (dotted line/asterisks) is compared to the revised 
fine-grained age distribution used here (dashed line/diamonds). 
The age bins of the spawned particles have widths of A log r = 0.5 
as in K10. 



Following the first TV-body simulation of the Anten- 
nae Galaxies by TT72, several groups have i mproved 
the pure stellar dynamical models of the system ([Barnes! 
fl98l IHibbardl l200l IRenaud et all 120081 [2009h . How- 
ever, there are still only a few ./V-body+hydrodynamical 
simulations that predict star formation histories of the 
Antennae and, at the same time, provide a reasonable 
match to the gross morphology and kinematics of the 
system. 

In the MBR93 simulation, the ISM is represented by 
discrete clouds that evolve by merging and fragmenta - 
tion. Star formation is modeled using a ISchmidtl <fl959h 
relation, cx p 18 , with an initial gas de pletion timescale 
of w 7Gyr (see lMihos et"allll99l fl99l for details). The 
spatial and mass resolutions for the baryonic component 
(stars) are x ICS = 200 pc and m res = 6.8 x 10 6 M Q . The 
Antennae w ere modeled a s two interacting identical disk 
galaxies (see iBarncs 1988|) on an elliptical orbit with ini- 
tial disk inclinations as in TT72. At the time of best 
match, after the first pericenter passage, the total SFR 
has increased by a factor of 6 with respect to the progen- 
itor disks. Most of the star formation is concentrated in 
the two nuclei, with very little in the overlap region be- 
tween them. Observations show, however, that most of 
the recent star and cluster formation takes place in the 
overlap region ie^.Stanford_et_aL| [1990; Mir abel et all 
Il998t IZhang et al.ll200lD . We have reproduced the star 
formation history from Figure 10 of MBR93 in our Fig- 
ure Q] (red dotted line). The intermittent nature of star 
formation in the MBR93 simulation, i.e. the rapid fluctu- 
ations in the overall SFR, is very likely a result of their 
particular discrete-cloud model of the ISM. To extend 
the simulated star formation histories to r > 10 9 yr, we 
assume a constant rate of star formation before the start 
of the simulations (i.e., the time at which the full N- 
body+hydrodynamical interactions are switched on) We 
adopt the same procedure for all the other simulations 
presented in the remainder of this section (see the con- 
stant horizontal lines in Figure [lj . 

Recently, K10 presented the first simulation of the 
Antennae that accounts for the extended star for- 
mation in the overlap region. For this study, the 
ISM was mo deled by smooth ed particle hydrodynam- 
ics (SPH; see [Monaghanl IT992L for a review). Star for- 
mation again followed a Schmidt relation, with cx 



,o 3 / 2 and a characteristic gas depletion timescale of 
to = 8.4 Gyr at the star formation threshold, n cr ; t = 
0.128 cm" 3 . The K10 simulation also includes stel- 
lar feedback (|Springel fc Hernquistl 120031) . The bary- 
onic spatial and mass resolutions are x res — 35 pc and 
irires = 7x 10 4 Mq. The progenitor galaxies were merged 
on a mildly elliptical orbit with disk orientations and 
a time of best match that resulted in better agreement 
with the observed large- and small-scale morphology and 
line-of-sight kinematics compared with previous Anten- 
nae simulations. In particular, the K10 simulation is the 
only published simulation that gives a reasonable spa- 
tial distribution of recent star formation activity, includ- 
ing the observed sta rburst in the overlap region (see also 
iKotarba et "all 1 20 101 ). 

K10 presented the stellar age distribution within the 
central 18 kpc of the simulation at the time of best match 
(see their Figure 4). This distribution was based on the 
stellar particles spawned from the gas particles accord- 
ing to the Schmidt relation. The time sampling for the 
particle spawning adopted by K10, however, is too coarse 
for the detailed analysis presented here. This can be seen 
in Figure [2j where we compare the total star formation 
rate (solid line) with the age distribution of spawned par- 
ticles (asterisks, 0.5 dex age bins). For ages older than 
a few x 10 7 yr, they are very similar. For younger ages, 
however, there are large fluctuations in the age distribu- 
tion of spawned particles. In particular, this distribution 
is too high in the first bin and too low in the second 
bin. We emphasize that these fluctuations are artifacts 
caused by the coarse time-sampling of the spawning pro- 
cess; hence they have only a numerical, not a physical 
origin. To remedy this inconsistency, we have rerun the 
K10 simulation with a finer time sampling for the particle 
spawning and show the corresponding age distribution 
(0.5 dex binning) in Figure [2] (open diamonds). Now the 
age distribution closely follows the SFR even for young 
ages (as it should) . In the following, we dispense entirely 
with the age distribution of spawned stellar particles and 
adopt the SFR computed directly from the Schmidt re- 
lation at each time step. 

For the present study, we analyze the simulation 
presented in K10 together with two new simulations 
with the same orbital and other model parameters, 
but with different amounts of feedback. In these 
simulations, stellar feedback heats and pressurizes the 
ISM, thereby suppressing and r egulating star forma- 
tion (jSpringel k, Hernquistl l2003f ). The energy input 
is controlled by a dimensionless parameter g cqs that 
ranges from (no feedback) to 1 (full feedback, see 
iSpringel et al.1 [20051 for details). In Figure [TJ we show 
the star formation histories for the simulations with full 
feedback (q cqs = 1, blue dashed), intermediate feedback 
(^oqs = 0.5, blue solid) from K10, and very weak feedback 
(<7cqs = 0.01, blue dot-dashed). The SFR in the simula- 
tion with the strongest feedback (q cqs = 1) increases only 
modestly after the first pericenter passage (vertical blue 
bar) and then again after the second pericenter passage, 
until it is 4 times higher at the time of best match than 
in the pre-interaction disk phase. Reducing the feedback 
(icqs = 0.5) results in a more efficient consumption of 
gas after the first close passage. At the time of best 
match, the SFR has increased only by a factor of 5. In 
the case of weak feedback (q cqs = 0.01), the gas is con- 
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Fig. 3. — Age distribution of star clusters derived from the star formation histories in the simulations shown in Figure [T] (without 
disruption ) . For comparison, we plot the observed age distribution of clusters with M > 2 X 10 5 Mq in the Antennae (filled dots) from 
I Fall et aLl ||2005I 1. The diagonal line represents the power-law model, dN/dr oc t 7 with the best-fitting exponent 7 = —1.01 ± 0.03. The 
simulated age distribution has been normalized to match the observed one at the midpoint of the youngest bin. The zeropoint of age is the 
time of best match in the simulations. The colored arrows indicate the starting times of the simulations. 



sumed more rapidly in an early interaction phase, with a 
peak in star formation about 150 Myr after the first peri- 
center passage. Thereafter, the SFR drops rapidly and 
then increases following the second pericenter passage to 
5 times the initial value at the time of best match. Most 
of the star formation then occurs in the overlap region 
rather than in the two nuclei, in reasonable agreement 
with observations (["Zhang et atll2001t IWang et all 120041 : 
iBrandl et al.ll2009t iKlaas et al.ll2010l ). 

TCB10 have recently presented a set of Antennae sim- 
ulations based on the adaptive mesh refinement (AMR) 
grid code RAMSES (iTevssierl f2002h with the o rbit of 
an e arlier iV-body simulation by iRenaud et all (|2008l 
2009). They model the ISM with a polytropic equa- 
tion of state. Stellar feedback is not included, but 
thermal support is added on small scales to avoid ar- 
tificial fragmentation. As in the simulations described 
above, star formation is modeled using a Schmidt re- 
lation in the form oc p 3 / 2 . TCB10 present low- 
resolution (x res = 96pc,m res = 10 6 Mq) and high- 
resolution (x res = 12pc,m res = 4x 10 4 Mq) simulations 
in which the star formation thresholds are scaled (by al- 
most two orders of magnitude higher in the latter case) 
to give similar SFRs in the progenitor disks. The gas 
depletion timescales at the star formation threshold are 
16 Gyr and 2 Gyr, respectively, for the low- and high- 
resolution simulations. The star formation history of the 
TCB10 low-resolution simulation is shown by the green 
dashed line in Figure [TJ As in the MBR93 simulation, 
the SFR increases moderately after the first pericenter 
passage (vertical green bar) , by a factor of 4 at the time 
of best match. In the high-resolution simulation (green 
dot-dashed), the increase in the SFR after first peri- 
center passage is more dramatic, about 20 times higher 
at best match than in the progenitor disks. The star 
formation history, however, parallels the low-resolution 



case with an almost constant offset. In the K10 simu- 
lation with q cqs = 0.01, a similar rise in the SFR after 
the first pericenter passage resulted from weak feedback. 
The high-resolution TCB10 simulation shows a more ex- 
tended distribution of star forming sites, but these are 
not as smoothly connected as in the K10 simulations or 
as in the real Antennae (see Figure 1 in TCB10). Fur- 
thermore, the TCB10 simulations do not reproduce the 
observed starburst in the overlap region. 

The star formation histories in the simulations assem- 
bled here show a great deal of diversity, due to different 
orbits (merger timescales, orientations, etc) and/or dif- 
ferent prescriptions for star formation and stellar feed- 
back. The differences in the SFRs can be as much as an 
order of magnitude at a given evolutionary phase (e.g., 
the first pericenter passage). In all cases, the SFR is en- 
hanced by the interaction relative to the quiescent disk 
phase, by factors of 4 — 20 at the time of best match. 
However, most of this variation occurs early during the 
interaction, shortly after the first pericenter passage. In 
contrast, for r < 2 x 10 8 yr, the SFR varies remarkably 
slowly, by factors of 3 or less. This near-constancy of the 
recent SFR plays an important role in our interpretation 
of the observed age distribution of the star clusters in 
the Antennae galaxies. 

4. COMPARISONS WITH OBSERVATIONS 

We now compare the simulations with observations. 
The smooth curves in Figure [3] show the same simulated 
SFRs as FigureQ] but now plotted as \og(dN/dr) against 
log(r/yr) with the zeropoint of age (r = 0) taken to be 
the time of best match for each simulation. The data 
points in Figure [3] show the observed age distribution of 
massive clusters (M > 2 x 10 5 Mq) in the Antennae from 
iFall et all ()2005f ). The statistical uncertainties (VN) in 
each bin are 0.05 - 0.06 dex, about half the radius of the 
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Fi g. 4. — Implied survival fraction of star clusters in the Antennae, / SU rv = {dN / dr) / (dN / dr)f ovui , plotted against age, with dN/dr 
from Fall et al. (2005) and (dN/dT){ orm from the simulations shown in Figure [T] (binned in the same way as the observations). The colored 
arrows indicate the starting times of the simulations. 



plotted symbols. The diagonal line in Figure [3J shows the 
power law that best fits the observations: dN/dr oc r 7 
with 7 = —1.01 ± 0.03. We have normalized the simu- 
lated and observed dN/ dr to a common value at the first 
data point. What is striking about Figure [3] is the di- 
vergence between the simulations and observations with 
increasing age, indicative of progressive disruption. At 
r = 10 s yr, for example, the observed age distribution 
has declined by a factor of 42, while the simulated forma- 
tion rates have declined by factors of only 1.3 - 2.5, even 
in the midst of a dramatic collision between the galaxies. 
This illustrates the main conclusion of this paper: both 
the observed age distribution and simulated formation 
rates of clusters in the interacting Antennae galaxies re- 
semble those in quiescent (non-interacting) galaxies. We 
demonstrate this conclusion more quantitatively in the 
remainder of this section. 

The age distribution shown in Figure UJ is based on 
UBVIHa images taken with the Wide Field Planetary 
Camera 2 on HST. The masses and ages of the clusters 
were derived by comparing their luminosities and colors 
with stellar population models while correcting for inter- 
stellar extinction. The resulting age distribution is not 
sensitive to the particular choice of models or extinction 
curve. The data points plotted in Figure [3] pertain to 
a mass- limited sample with M > 2 x 10 5 M Q , which is 
essentially complete for ages up to r m 10 9 yr. The age 
distributions for samples limited at lower masses have 
the same power-law shape but do not extend to such 
large ages. In constructing the age distribution, rela- 
tively large bins were chosen (Alogr = 0.8) in order to 
smooth over small-scale features caused by the system- 
atic errors that arise whenever masses and ages are esti- 
mated from multiband photometry. The most insidious 
of these is the artificial gap atrw (1—2) x 10 7 yr resulting 
from stochastic variations in the colors of clusters caused 
by rapid evolution of a few bright red supergiant (RSG) 



stars. This so-called RSG gap, and other, less prominent 
features, affect every age distribution in this field — not 
just the one for the Antennae. The only practical way to 
minimize the corresponding systematic errors in the age 
distribution is to choose relatively large bins and to avoid 
centering any of them on the RSG gap. In the future, 
it may be possible to reduce the prominence of the RSG 
gap by augmenting the UBVIKa photometry with addi- 
tional infrared photometry in several bands and/or spec- 
tra of larg e samples of clust ers. The interested reader is 
referred to lFall et~a l. (2005) for a more complete discus- 
sion of the age distribution shown in Figure [3l including 
the many tests that were performed to determine its ac- 
curacy and robustness. Confirm ation and further analy- 
sis ca n be found in the pap ers by IWhitmore et al.l (|2007l 
l20inh and lFall et all IpOOfl) . 

The common normalization of dN/ dr at small r in Fig- 
ure |3] for both simulations and observations is required 
by our assumption that the formation rates of stars and 
clusters are linked together by Equation ([2]). In principle, 
we should equate dN/dr and (dN/dr)f olm at r = 0, the 
only age at which we can be certain that no clusters have 
yet been disrupted, that they are all intact and observ- 
able, even if some of them are not gravitationally bound 
and will eventually disperse. This procedure implicitly 
fixes the constant of proportionality c in Equation ([2j be- 
tween the formation rates of stars and observed clusters 
(in this case, those with M > 2 x 10 5 M ). The dispersal 
time for unb ound clusters is At ~ 10 yr (~ 10 cross- 
ing ti mes: see lFall et al.ll2005l and lBaumgardt fc Kroupal 
l2007t ). Thus, by matching dN/dr and (dN/dr) {orm at 
r = 2.5 x 10 6 yr rather than r = 0, we likely make 
only a small error, of order unity or much less. It 
would be a mistake to equate dN/dr and (dN/dr){ OTU1 
at r >> 10 7 yr because in general this would cause a 
mismatch at young ages, the only part of the distribu- 
tion for which disruption is guaranteed to be negligible. 
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Fig. 5. — Predicted age distribution of star clusters based on the simulated formation rates shown in Figure[T]and the power-law model for 
the survival fract ion, / sur v oc t s , w ith the best-fitting exponents listed in Table [T] For comparison, we show the observed age distribution 
(filled dots) from Fall ct al. (2005). In the upper panel, the formation rates are unbinncd; in the lower panel, they are binned in the same 
way as in the observations. The colored arrows indicate the starting times of the simulations. 



We now explore several quantitative measures of the 
progressive disruption of clusters indicated by Figure OS 
First, we divide the observed age distribution by the 
simulated formation rates (after binning in the same 
way) to obtain the implied survival fraction, / surv = 
(dN/ dr) / (dN / dr) f olm . This is plotted as a function of 
age in Figure |U Evidently, / surv declines monotoni- 
cally up to t = 10 s yr for all simulations and up to 
r = 6 x 10 s yr for all but one simulation. The statistical 
significance of the decline in / surv is 17-20 a from the 
first bin to the second, 5-9 a from the second bin to the 



third, and 8-16 a from the third bin to the fourth (again 
with one exception). Thus, a constant / SU rv is definitely 
ruled out. A corollary of this result is that disruption 
continues at least up to r ~ 10 s yr and possibly up to 
t ~ 10 9 yr. From Figure 0] it also appears that the 
survival fraction declines roughly as a power law of age. 

We can make a more direct test of the power-law model 
for the survival fraction, / SU1V oc t s , by combining it with 
the formation rate (c?A/dr)f orm as in Equations ([1]) and 
((2]) ■ We then obtain the best- fit values of S by minimizing 
X 2 between the predicted and observed age distributions 
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TABLE 1 

Best-fit exponent <5 in the power-law model of the survival fraction. 



Model 


K10 q = 1.0 


q = 0.5 


q = 0.01 


TCB10 low-res 


hi-res 


MBR1993 


<5 


-0.86 


-0.89 


-0.88 


-0.80 


-0.62 


-0.77 


Errors 


0.01 


0.01 


0.02 


0.01 


0.02 


0.01 



dN/dr. The results of this comparison are listed in Table 
[T] and shown in Figure [5] In the upper panel of the fig- 
ure, we plot the predicted dN/dr without binning, while 
in the lower panel, we plot it with the same binning as 
the observations. The best-fit exponents lie in the range 
—0.9 < (5 < —0.6 (with minor differences among the 
simulations). These are similar to, but slightly larger, 
than the exponent 7 = —1.0 of the observed age dis- 
tribution, because the simulated formation rates decline 
gradually with age (at least for r < 2 x 10 s yr). We 
do not expect perfect agreement between any of the pre- 
dicted and observed age distributions, because, as noted 
in Section[3l there is substantial diversity among the sim- 
ulated formation histories themselves. Nevertheless, as 
Figure [5] shows, the predicted age distributions for all six 
simulations follow the same general, power-law trend as 
the observed age distribution. This is additional support 
for our conclusion that clusters in the Antennae are dis- 
rupted progressively over an extended period and in a 
manner similar to that in more quiescent galaxies such 
as the Milky Way and the Magellanic Clouds. Is is also 
gratifying to note that the simulation that best matches 
the age distribution (the K10 Qeqs = 0.01 simulation, 
dot-dashed blue line in Figured]) also reproduces the ex- 
t ended star formation in the overlap region. 

Bast ian et al.l (|2009T ) reached a different conclusion: 
that the disruption of clusters ceases at r ~ 10 7 yr in 
the Antennae. Th eir analysis di ffers from ours in sev- 
eral respects. (1) iBastian et al.l base their claim only 
on the MBR93 simulation. This simulation, however, 
fails to reproduce the observed spatial extent of recent 
sta r formation in the overlap region of the Antennae. 
(2) IBastian et al.l equate the observed age distribution 
dN/dr and the simulated formation rate (dN / dr)i mm at 
relatively old ages, r ~ 10 s yr. As we have emphasized 
above, dN/dr and (dN/dr)f orul must be matched near 
r « in order to connect the fo rmation rates of clus- 
ters and stars. (3) IBastian et al.l have adopted an age 
distribution with a dip at r = (1 — 3) x 10 7 yr and a sec- 
ondar y peak at r = (3—10) x 10 7 yr from lWhitmore et al.l 
(|2007[ ). This is based on the same observations as the 
age distribution adopted here, but it has narrower bins 
(Alogr = 0.5), one of which, unfortunately, is centered 
right on the RSG gap 4 . This dip-peak structure in the 
age distribution is not reproducible by any of the simu- 
lations considered here and we therefore regard it as an 
unphysical artifact associated with the RSG gap. 

5. CONCLUSIONS 

The main conclusion of this paper is that the star clus- 
ters in the interacting Antennae galaxies are disrupted 



in much the same way as those in other galaxies. In 
most if not all star-forming galaxies, including the An- 
tennae, the observed age distribution of clusters can be 
approximated by a power law, dN / dr oc r 7 , with an ex- 
ponen t in the range -1-0 < 7 < -0.7 fsee lChandar et all 
l2010bl and Section [2] here) . In general, this must reflect 
the combined formation and disruption histories of the 
clusters. However, variations in the formation rate are 
expected to be a minor influence because dN/dr is so 
similar in different galaxies and declines by such a large 
factor, typically ~ 10 2 , over a relatively small range of 
a S e i r ^ 10 8 -10 9 yr (i.e., less than 10% of the lifetime 
of the galaxies). Indeed, in several well-studied galaxies 
(the Milky Way and the Magellanic Clouds), the SFR 
is known from observations to have been nearly con- 
stant for the past ~ 10 9 yr, compelling evidence that 
the d ecline in dN/ dr is mainly a consequence of disru 
tion dBinnev fc Merrifieldlll998tlHarris fc Zaritskvll200 
20091 IChandar et al.ll2010al) . 



4 In fact, Whitmorc ct al. (2007) present two age distributions 
(in th eir Figures 4 and 15). The first one, adopted bv lBastian et af] 
(2009), was intended mainly for comparison with Monte Carlo sim- 
ulations. When we fit a power law to this age distribution, we 
ob tain 7 = —0.98 ± 0.14. The second age distribution presented 
by IWh itmorc et al. is essentially the same as the one presented by 
IFall et al.l (120051 ). 



The interpretation is less straightforward for the An- 
tennae galaxies, since they are too far away to determine 
their star formation history from observations. Further- 
more, it is natural to wonder whether the interaction 
could trigger enough recent formatio n to explain the 
shape of dN/dr without disruption. IFall et all ()2005f ) 
argued that the formation rate would vary by factors of 
a few on the orbital timescale (~ 10 8 yr), too gradu- 
ally to_accourrt_Jor mo&t of the_de£hne in_d/V/dr (see 
also lWhitmore et al.ll2007l and lFall et al.ll2009D . The re- 
sults presented here support this suggestion. We have as- 
sembled the star formation histories in all the available 
-/V-body+hydrodynamical simulations of the Antennae. 
These are based on different numerical methods, differ- 
ent orbits, and different prescriptions for star formation 
and stellar feedback. The treatment of small-scale pro- 
cesses is still approximate at best, due to the low reso- 
lution in the simulations compared to molecular clouds 
and clumps in the real ISM. The star formation rates 
differ greatly among the simulations, both in absolute 
level and in time dependence. Nevertheless, we find that 
they all vary slowly, by factors of only 1.3 — 2.5 in the 
past 10 8 yr. When we combine the formation rates in 
the simulations with a power-law model for the survival 
fraction of clusters, / SU rv oc t 5 , we find good agreement 
with the observed age distribution over the range 10 6 yr 
^ r < 10 9 yr for -0.9 < 5 < -0.6. The similarity be- 
tween S and 7 indicates that dN/dr is shaped mainly by 
the disruption of clusters rather than variations in their 
formation rate. 

The only caveat to this conclusion stems from our as- 
sumption that the formation rates of clusters and stars 
are proportional to each other, i.e. {dN/dr)f orrn — 
c ■ dN*/dr with c = constant (cf. Equation This is 
certainly true if most stars form in clusters, a hypoth- 
esis consistent w ith Ha observations of the Antennae 
(jFall et al.M2005h . However, even if we were to abandon 
this assumption entirely, and allow c to vary arbitrarily 
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with r, a slightly modified version of the analysis pre- 
sented above would then lead to the alternative result 
c(t) • / surv (r) oc t s with -0.9 < 5 < -0.6. We would 
then be left with the problem of explaining why the prod- 
uct c • / sur v but not / surv in the Antennae happens to be 
the same as / SU rv alone in other galaxies. The simplest 
interpretation — the one advocated here — is that the 
formation rates of clusters and stars do track each other 
and that the survival fractions and hence the disruption 
histories are similar in different galaxies, including the 



Antennae. 
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